hep-ph/9612444l Revised EHU-96-12 
June 1997 



Prediction on CP Violation in Long Baseline 
Neutrino Oscillation Experiments 



Morimitsu TANIMOTO Q 

Science Education Laboratory, Ehime University, 790-77 Matsuyama, JAPAN 



ABSTRACT 

We predict CP violation in the long baseline accelerator experiments taking into 
consideration the recent LSND data and the atmospheric neutrino data. The estimated 
upper bound of CP violation is 0.006, which may be observable in the long baseline 
accelerator experiments. It is found that the upper bound increases to 0.01 if the LSND 
data is excluded. The matter effect, which is not CP invariant, is found to be very 
small in the case we consider. 



E-mail address: tanimoto@edserv.ed.ehime-u.ac.jp 



1 Introduction 



The origin of CP violation is still an unsolved problem in particle physics. In the 
quark sector, CP violation has been intensively studied in the KM standard model. 
UU For the lepton sector, CP violation is also expected, unless neutrinos are massless. 
If neutrinos are Majorana particles, one finds three CP violating phases. [0 However, 
the effect of extra Majorana phases is suppressed by the factor {m v / E) 2 . || Therefore, 
CP violation in the neutrino flavor oscillations relates directly to the CP violating 
phase in the mixing matrix for massive neutrinos. 

Unfortunately, this CP violating effect is suppressed in the short baseline neutrino 
oscillation experiments if the neurinos have a hierarchical mass spectrum. However, 
suppression is avoidable || |||] in the long baseline accelerator experiments which are 
expected to operate in the near future. 0, || Observation of the CP violating effect 
may be possible in those experiments. 

The recent indications of a deficit in the flux of the atmospheric neutrinos |§ 
has renewed interest in using accelerator neutrinos to perform long baseline neutrino 
oscillation experiments. Many possibilities of experiments have been discussed. [0, [8| 
The first experiment will begin in KEK-SuperKamiokande. || A few authors have 
expressed interest in CP violation in those experiments. || |6| We also discussed CP 
violation and estimated its upper bound. However, we did not take into account all 



available data, for example, the recent LSND data in our previous estimate. By 
using the LSND data, we can obtain unambiguous results for CP violation. Thus, 
the purpose of this paper is to predict CP violation in those accelerator experiments 



taking into consideration both the recent LSND data [12] and atmospheric neutrino 



data. [11] 

On the other hand, the background matter effect, which is not CP invariant, is 
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a significant obstacle in the attempt to observe CP violation originating from the 
phase of the neutrino mixing matrix. We present careful analyses for CP violation 
in the background matter. Therefore, we give in this paper more definite predictions 
compared with the previous ones. M 

In §2, we discuss the neutrino mass pattern, in which the recent LSND data |12| and 



atmospheric neutrino data|TT[ are consistent with the predicted transition probabilities. 
In §3, the general framework of CP violation and the matter effect are discussed. We 
predict the magnitude of CP violation including the matter effect in §4. Section 5 is 
devoted to a conclusion. 

2 Neutrino mass hierarchy and LSND data 

One expects to extract the neutrino masses and flavor mixing matrix from the data of 
neutrino flavor oscillations. Although recent short baseline experiments have not yet 
confirmed the evidence of neutrino oscillation, they suggest some patterns of neutrino 
masses and mixings. For example, the data given by the LSND experiment |12| indicate 



at least one large mass difference, such as Am 2 ~ 0(1 eV 2 ) with sin 2 2#lsnd — 10 3 . 



In the near future, we can expect the data from the KARMEN experiment [13] which is 
also searching for the — > z/ e (z/ M — > V e ) oscillation as well as LSND, and the data from 
the CHORUS and NOMAD experiments, [13] which are looking for the — > v T oscilla- 



tion. The most powerful reactor experiments searching for the neutrino oscillation are 
those of Bugey |1| and Krasnoyarsk |16[ at present. They provide excluded regions in 



the (sin 2 29, Am 2 ) parameter space by non-observation of the neutrino oscillation. The 



first long baseline reactor experiment CHOOZ [17 will also give us severe constraints. 



These experiments may determine patterns of neutrino masses and mixings. For the 



present, constraints are also given by the disappearance experiments in CDHS |L8 and 
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CCFR, pi and appearance experiments E531, |2(J E776 Q and new CCFR. pf By 



use of those experimental bounds, some patterns are studied for neutrino masses and 
mixings. p| |7| 

In this paper, we take account of the LSND data seriously. Okada and Yasuda 
presented clearly [^8| that the Bugey reactor experiment |15[ and the CDHS |18| ac- 



celerator experiment exclude the mass hierarchical pattern m 3 3> m 2 3> m^or ~ mi) 
if the LSND data are used. Therefore, the mass pattern 772,3 — m 2 ^> mi is favored. 

If the excess of the electron events at LSND are due to the neutrino oscillation, 
the v ^ — > v e oscillation transition probability is equal to 0.311qjo 

±0.05 % [0. The 

plot of the LSND favored region of Am 2 vs sin 2 29 allows the mass region to be larger 
than Am 2 ~ 0.04 eV 2 . However, the constraint by Bugey [15| has excluded a mass 



region lower than Am 2 = 0.27 eV 2 . On the other hand, the constraint by E776 |21] 
has excluded a region larger than Am 2 = 2.3 eV 2 . Thus, the mass squared difference 
Amfi is obtained in the range 

Am 2 x ~ 0.27 eV 2 ~ 2.3 eV 2 . (1) 

There are only two hierarchical mass difference scales Am 2 in the three-flavor mixing 
scheme. If the highest neutrino mass scale is taken to be O(leV), which is appropri- 
ate for cosmological hot dark matter (HDM), |29| the other mass scale is either the 
atmospheric neutrino mass scale Am 2 ~ 10~ 2 eV 2 |TTJ or the solar neutrino mass 
scale Am 2 ~ 10~ 5 ~ 10~ 6 eV 2 . |0| Since the long baseline experiments correspond 



to the atmospheric neutrino mass scale, we take Am 2 ~ 10 2 eV 2 as the lower mass 



scale. Then, one should introduce the sterile neutrino [31] if one wishs to solve the 
solar neutrino problem from the standpoint of a neutrino oscillation. 

Since we fixed the neutrino mass scales, we can now discuss the pattern of the 3x3 
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neutrino mixing matrix. In general, there are three allowed regions of the mixings, 
which are derived from the reactor and accelerator disappearance experiments. 

The Bugey reactor experiment JOJ and CDHS [[L8[] and CCFR |T9| accelerator 
experiments give bounds for the neutrino mixing parameters at the fixed value of 
Amjj for three cases, as seen in Appendix A: 



(A)\U el \~l, |C/ M i|<l, (B)\U el \<l, |M<1, {C)\Uei\<l, \U^\~1. 

(2) 

The cases (A) and (C) are consistent with the LSND data because |Z7 e i^ti| could 
be O(10~ 3 ), as seen in Eq. (A- 6). [?§ On the other hand, the cases (A) 
and (B) are consistent with the atmospheric neutrino data. The following transition 
probability formulae for the atmospheric neutrino may be helpful to understand the 
above situation: 



- "r) = -WvUrtU&Urt sin 2 (^|^) + 2U^U^ , 
PW ^ »*) = -W e2 U„ 2 U e3 U, 3 sin 2 (^J^) + . (3) 

We calculate the CP violation effect in the following pattern of the mixing matrix 
(case (A)), which is only consistent with both LSND and atmospheric data: 



1 


^3 






£2 


u r2 



\U\* I 6i U, 2 U, 3 I , (4) 
where £j (i — 1 ~ 4) are small numbers. 

3 CP violation and matter effect 



The direct measurements of CP violation originating from the phase of the neutrino 
mixing matrix are the differences of the transition probabilities between CP-conjugate 
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channels (see Appendix B): 

AP = P{v^ v e ) - P(^ -»• v e ) = P{v^v T )-P{V^V T ) 

= P(v e -> 17 r ) - P{y e -> i/ T ) = 4J£ P (sinPi 2 + sinP 2 3 + sinP 3 i) , (5) 



where 



,2 ^ 



Aj = AmJ— , (6) 



and |2 



J£p = Im^^l^gl^-jC/,^) = Si2S23Si 3 Ci2C23Ci 3 sin0 . (7) 

In the short baseline experiments, we obtain P 31 ~ — D 12 = 0(1) if Amjj = 
0(1 eV 2 ). The factor (sinPi 2 -t-sinP 2 3 + sinP 3 i) is suppressed because the two largest 
terms almost cancel due to their opposite signs. Another term is still small. Therefore, 
the CP violating quantity in Eq. (|5|) is significantly reduced. So one has no chance to 
observe CP violation for the present in the short baseline neutrino oscillation experi- 
ments. However, the situation is different in the long baseline accelerator experiments. 
The oscillatory terms sin P 12 and sin P 31 can be replaced by the average value since 
the magnitude of P 3 i(— — Di 2 ) is 10 3 ~ 10 4 . Then CP violation is dominated by the 
sinL>23 term, which is 0(1). Thus the CP violating quantity AP is not suppressed 
unless Jq P is very small. 

In order to predict the magnitude of CP violation, one should estimate Jq P - In 
case (A), which is consistent with the LSND data and the atmospheric neutrino data, 
we should take C12, Ci 3 ~ 1 and C2 3 — l/v2- The LSND data determine the value of 

Ml (R\ 



c 12 c 13 | c 23 s 12 + S23<Sl3Ci2e 



as seen in Eq. (A- 6), if Am^ is fixed. In the future, the observation of the process 
v e — > v T or — > v T will give a constraint for 

ci2Ci 3 |s 23 si2 - c 23 si 3 ci 2 e^| . (9) 



On the other hand, the disappearance experiment in the work by Bugey[jT5| gives 
a bound for C12C13 which is not so confining. For example, we obtain: |24 



|ci 2 ci 3 | > 0.984 at Am 2 31 = 2 eV 2 . (10) 

Using these constraints, we find the allowed parameter regions of su, S13 and <fi. 

We show the allowed curves in the S12 — s 13 plane for the fixed values of <fi and 
s 23 = l/v2 in Fig. 1, in which we have taken the typical LSND data: Am^ = 2eV 2 
and sin 2 2#lsnd — 2.6 x 10~ 3 . It is easily understood that larger s 12 and S13 are allowed 
as cf) increases to 180°. As seen in Eq. (^J), the larger S12 and S13 near 180° are expected 
to lead to a larger transition probability of v e — > v T or — > v T . The predicted 
maximal probabilities at CHORUS and NOMAD are P{y e — > u T ) ~ 3.3 x 10~ 4 and 
P{ v n — * v t) — 2.6 x 10 -8 , which are still below the sensitivity of those experiments. 

It may be useful to comment that the atmospheric neutrino anomaly is attributed 
to the — ► v T oscillation due to S23 = l/y/2 in this case. 

Now, we discuss the matter effect in the case of m 3 ~ m 2 ^> m\. The general 
discussion of the matter effect in the long baseline experiments was given by Kuo and 



Pantaleone. |3^] The T violation effects in the earth were also studied numerically by 
Krastev and Petcov. Their results indicate that there exist relatively large regions 
of masses and mixings for which the T violation effect can be considerably enhanced 
due to the presence of matter when neutrinos cross the earth along its diameter. 

Although the distance travelled by neutrinos is less than 1000 Km in the long 
baseline experiments, those data include the background matter effect which is not 
CP invariant. The matter effect should be carefully analyzed since the effect depends 
strongly on the mass hierarchy, mixings, and the incident energy of the neutrino, as 
shown in previous works. [[33|, |3^ 



The analytic formulation of the matter effect is given for the three generation case in 



Appendix C. [|35| The matter effect in the long baseline accelerater experiments is rather 
easily calculated by assuming a constant electron density, which is n e = 2.4 mol/cm 3 . 
However, the mixing angles and the CP phase, in matter depend in a quite complicated 
way on the vacuum mass-squared differences, vacuum mixing angles and vacuum CP 
phase as shown in Appendix C. 

We give qualitative discussion about the matter effect by using the approximate 
formula of the effective mixing angle 9 m i2, which is given in terms of vacuum mixings 



as 



33 



. oa M sin 29 12 

sin20 m i2 = = = , (11 



'(A cos 2 13 - Am 2 ,! cos 26 12 ) 2 + Am^ sin 2 26 12 
to zeroth order in Asm28i 3 . In the present case of m 3 ~ m 2 = 0(1 eV) ^> m 1; we 

have the relation Am^ > A ~ 5 x 1(T 4 eV 2 (A = 2\/2G F n e E). Then, one can easily 

find that the the matter effect A contributes slightly to sin# m i 2 in Eq. (|TT]) as long 

as 6 12 *C 1. The matter effect on the angles # m23 and # m i 3 are also expected to be 

very small. Thus the matter effect is expected to be small in the case we consider. 

Numerical studies are discussed in the next section. 

4 CP violation in long baseline accelerator experi- 
ments 

Long baseline accelerator experiments are planed to operate in the near future. J?| |[ 
The most likely possiblities are in KEK-SuperKamiokande (250 Km), CERN-Gran 
Sasso (730 Km) and Fermilab- Soudan 2 (730 Km) experiments (MINOS). The aver- 
age energy of the beams are approximately 1 GeV, 6 GeV and 10 GeV at KEK- 
PS (12 GeV), CERN-SPS (80 GeV) and the Fermilab proton accelerator (120 GeV), 
respectively. 

We can estimate the factor sinZ?i 2 + sin/} 2 3 + sinD3i in those experiments. || 
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Although estimated values depend on Ara^ and E, these could be close to 1 for 
Anijj = 1 ~ 10 eV 2 . In order to predict the magnitude of CP violation at the KEK- 
SuperKamiokande (250 Km) experiment, we fix L = 250 Km and take E — 1.3 GeV 
with Am| = 0.01 eV 2 . This experiment will be a disappearance experiment, since such 
neutrinos will not produce tau leptons. But there is a manifestation of CP violation 
even in the case of disappearance if final states are summed over. 

Now we can calculate AP(z/ At — > v T ) = P(z^ M — > v T ) — Piv^ — > v T ) and APfV^ — > 
u e ) = P{ffj, — > v e ) ~ Ptyfj. ~~ > v e ), which are the direct measurement of CP violation 
in the long baseline experiments. In the present case, Piy^ — > v T ) is large (~ 0.4) but 
P( v n ~ * v e) is very small (~ 0.01). By use of the constraints in Fig. 1, we predict 
AP(z/ M — > u T ) and AP(z/ M — > u e ) in the long baseline experiment with L = 250 Km and 
E = 1.3 GeV. We show AP in the matter versus the vacuum mixing s\2 with <p = 159° 
in Fig. 2, in which the predicted AP(z/ M — > u T ) in the vacuum is also presented by the 
dotted curve. The predicted AP(z/ M — > v e ) in the vacuum cannot be distiguished from 
that in matter. It is found that the matter effect is very small. As <fi changes from 
90° to 180°, larger s\2 and S13 are allowed, but sin0 decreses. So the maximal Jq P 
is realized between 90° to 180°. Actually, the largest value of CP violation is given 
at (f) = 159°. The predicted largest value is 0.006, which may be obserbable in the 
KEK-SuperKamiokande experiment in the future. 

It is noticed that the CP non-invariant quantity |AP(^ — > u T )\ is different from 
|AP(z/ At — > u e ) \ in matter. On the other hand, these have the same magnitudes for CP 
violation originating from the phase in the neutrino mixing matrix, as seen in Eq. (||). 

In our predictions, we have used the LSND data. However, it may asserted that 
the LSND result has not yet been confirmed. One has to wait for the data from the 
KARMEN experiment which is also searching for the — > z/ e (l7 M — ► V e ) oscillation 
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as well as LSND. In order to clarify the role of the constraint from the LSND data, 
we discuss the case excluding the LSND data. As long as the highest neutrino mass 
scale is taken to be 0(1 eV), which is appropriate for the cosmological HDM, and the 
atmospheric neutrino anomaly is due to the large z/ M — > v T oscillation, the mass pattern 
777.3 — m 2 ^ m i is also favored, even if the LSND data is excluded (here, the hierarchy 
of masses is assumed). Then, the most significant constraint is given by the E776 data. 
The mass squared difference Am^ is allowed in the region larger than 2.3 eV 2 , which is 
the upper bound given by LSND. In the range Am^ = 2 ~ 10 eV 2 , we find the largest 
bound of sin 29 = 3 x 10~ 3 at Am§! = 8 eV 2 , which is not allowed by the LSND data. 
By use of this bound, we show the curves of the upper bounds in the S12 — s 13 plane, 
fixing S23 = l/v2 in Fig. 3, where the larger of s±2 and S13 are allowed. The largest 
value of the CP violation is given at <fi = 163°. Predicted curves of the upper bounds 
in the matter are shown in Fig. 4. The predicted largest value is 0.01 if the LSND data 
is excluded. Thus, the LSND data provides a somewhat severe constraint. 

5 Conclusions 

We have studied the direct measurements of CP violation originating from the phase of 
the neutrino mixing matrix in the long baseline neutrino oscillation experiments. We 
can predict CP violation in these accelerator experiments taking into consideration 



both recent LSND data |L2] and atmospheric neutrino data. |TTJ These data favor the 
mass hierarchy 7773 ~ 7772 ^> 777 1 with 777.3 = 0(1 eV) and the flavor mixings \U e \\ ~ 1 
and I C/^i I <C 1 with jt/^l — l/v^- The estimated upper bound of CP violation is 
0.006, which may be observable in the long baseline accelerator experiments in the 
future. The matter effect, which is not CP invariant, is found to be very small. 

In order to clarify the role of the constraint from the LSND data, the case excluding 
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the LSND data has been discussed. Then, the predicted upper bound is 0.01. It is 
noted that the LSND data provides the somewhat tight constraint. 
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A Experimental Constraints 

We show constraints for the neutrino mixings from the reactor and accelerator disap- 
pearance experiments. Since no indications in favor of neutrino oscillations were found 
in these experiments, we only obtain the allowed regions in (U^, Am^) parameter 
space, where Am| = mf — m| is defined. The Bugey reactor experiment |I5| and 



CDHS and CCFR |19| accelerator experiments give bounds for the neutrino mix- 
ing parameters at the fixed value of Am^. We follow the analyses given by Bilenky et 
al.El 



Since the CP violating effect can be neglected in those short baseline experiments, 
we use the following formula without CP violation for the probability in the disap- 
pearance experiments: 

P(u a - u a ) = 1 - 4|f/ Ql | 2 (l - \U al \ 2 ) sin 2 (^f^) , (12) 

where the mass hierarchy Am^ ~ Am^ ^> Am^ is assumed. The mixing parameters 



can be expressed in terms of the oscillation probabilities as p4 

1 



U al \ 2 = -(l±^l-B UaUa ) , (13) 
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with 



B„ 



{l-PK-, Q )}sin- 2 (^|^; 



(14) 



where a = e or /i. Therefore the parameters U^i at the fixed value of Am^ should 
satisfy one of the following inequalities: 



\U al \ 2 > + B UaK 



,(+) 



or \U al \ 2 <-(l- Jl-B UaVa ) = aP 



(15) 



The negative results of Bugey, CDHS and CCFR experiments have given the values of 
and af\ @] For example, one obtains = 0.984, = 0.016, aj+> = 0.985 
and = 0.015 for Am^ = 2 eV 2 . 

It is noted from Eq. ([15]) there are three allowed regions of \U e i\ 2 and l^i] 2 as 
follows: 



(A) |t/el| 2 >4 + ) 

(5) \U el \ 2 <a^ 
(C) \U el \ 2 <a^ 



Ki\ 2 <ap, 
\U^\ 2 <a^, 



(16) 



\U, 



In addition to these constraints, we should take account of the constraints of E531 |20 
and E776 |2T| experimental data by using the transition probabilities 

.Am^I, 



P(i/„->i/ t ) ~4|^ 1 | 2 |f/ Tl | 2 sin 



AE ' 
.Am^I, 

■ AE ' 



(17) 



These constraints often become severer than the ones of the disappearance experiments. 



B CP Violation 



The amplitude of the v a —>■ v@ transition with the neutrino energy E after traversing 
the distance L can be written as 



A{y a -> vp) = e 



-iEL 



k=2 



exp 



2E 



1 
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where the U a i denote the elements of the 3x3 neutrino flavor mixing matrix, in 
which a and i refer to the flavor eigenstate and the mass eigenstate, respectively. The 
amplitude A{v a — > Up) is given by replacing U with U* on the right-hand side of 
Eq. fllTf). The direct measurements of CP violation originating from the phase of 
the neutrino mixing matrix are the differences of the transition probabilities between 
CP-conjugate channels: [Q 

AP EE P{V„ -> V e ) - P(Z/ M -> U e ) = P( Vfl -> Z/ T ) - P(I7 M -> F T ) 

= P(F e -> 17 T ) - P{y e -> i/ T ) = 4J£ P (sinP 12 + sinP 23 + sinP 31 ) , (19) 



where 



A, = Am J ^ , (20) 



and J^p is defined for the rephasing invariant quantity of CP violation in the neu- 
trino mixing matrix as well as that in the quark sector. 0] In terms of the standard 



parametrization of the mixing matrix, [EC 



C13C12 C13S12 si 3 e l<i, \ 

-C23S12 - S 2 3Sl3Cl 2 e^ C23C12 - S 2 3Sl3Si2e 4<?i S23C13 , (21) 
\ S23S12 — — S23C12 — C 2 3Si3Si 2 e^ C 2 3Ci3 / 

where = sin 8^ and Qj ee cos^- are vacuum mixings, and <fi is the CP violating 
phase, we have 

Jcp = lm(Uf, 3 U* 3 U* 2 U T 2) = si2S23Si 3 ci2C23Ci 3 sin0 . (22) 
The oscillatory terms are periodic in L/E, and P12 + P23 + P31 — is satisfied. 

C Matter Effect 

The effective mass squared in matter for the neutrino energy E in the weak basis 
35| is 



13 





{ m\ 





\ 




(A 







Ml = U m 










\ul + 








!) 




V o 





m\) 




vO 








(23) 

where A = 2\f2Gpn e E. For antineutrinos, the effective mass squared is given by 
replacing A — > —A and ?7 m — > L^. The effective mixing matrix U m is 



(24) 



— c m23 s ml2 — S m23 S m i 3 C m i 2 e l< ^ m C m 23 c ml2 — -S m 23 s ml3 s ml2e^ m S m 23 c ml3 
\ •S m 23S ?T i.l2 ~ Cm23Sml 3 Cml2e 4<7im — S m23 C m \ 2 — C m 23'Sml3'Sml2e i ^ m C m23 C m \ 3 / 

cos #mij are effective mixings, and m is the effective 



where s mij = sin mij and c mij 



phase in matter. 

If we use the constant electron density n e = 2.4 mol/cm 3 , the effective mixing angles 
and phase are given in terms of vacuum mixings and the effective neutrino masses Mi, 
M 2 and M 3 , which are eigenvalues in Eq. (f23|), as follows: p5| 



'ml2 



S 



where 



and 



■(M 2 — aM 2 + / 9)AM| 1 



AM| 2 (M 1 4 - acM? + (3) - AM| 1 (M 2 4 - «M 2 2 + /?) 
M 3 4 - aMf + p 



m13 AMlAMi 2 ' 

2 E 2 ^ + F 2 c 2 3 + 2EFc 23 s 23 cos0 

s m23 - ^ 2 + p2 ' ^ 

e - #m = (£ 2 e-* ~ F 2 e^) g2 3C 23 + £F(c 2 3 - s 2 3 ) 



(£ 2 S 2 3 + F 2 C 2 3 + 2EFc 23 S 23 cos0)(£ 2 s 2 3 + F 2 c 2 3 - 2EFc 23 s 23 COS0) 

22, 2/2 2 i 2 \ i 2/ 2 2 , 2 \ 

a = m 3 c 13 + m 2 (c 12 c 13 + s 13 ) + m^s^c^ + s 13 ) , 

,0 22/22, 2 2 \ i 222 

P = m 3 C l 3 ( m 2 C 12 + m l S 12) + ™2 m l S 13 > 

E = [Am^M 2 - m 2 2 ) - Am 2 21 (M% - m^s\ 2 ]c 13 s 13 , 

F = (Ml - ml)Am 2 21 c 12 s 12 c 13 , (26) 



Ml ~ i [m 2 + m\ + Ac 2 3 - \J(Ac 2 13 - Am 2 2l cos26 12 ) 2 + (Am 2 21 sin2^ 2 ) 2 ] 
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3 



ml + m 2 2 + Ac 2 13 + J (Acj 3 - Am^ cos 26» 12 ) 2 + (Am^ sin 29 12 ) 



Ml ~ m^ + A4. (27) 
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Figure 1: 

Allowed curves in the Si 2 — s i3 plane, fixing s 2 3 = 1/a/2, in the case Am^ = 2 eV 2 
and sin 2 2# L snd = 2.6 x 1CT 3 . The thick-solid curve, thin-solid curve, short-dashed 
curve, dash-dotted curve, dotted curve and long-dashed curve corrspond to 0=160°, 
150°, 120°, 90°, 60° and 30°, respectively. 
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Figure 2: 

Predicted AP(^ — > v T ) (solid curve) and AP(i/ /1 — > z/ e ) (dashed curve) in the matter 
versus s 12 with = 159°. AP(f /1 — > z/ T ) in the vacuum are also presented by the dotted 
curve. Here, E = 1.3 GeV and L = 250 Km with Am 2 32 = 1(T 2 eV 2 are used. 
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Figure 3: 

Upper bounds in the Si 2 — S13 plane, fixing s 2 3 = l/\/2, in the case Am^ = 8 eV 2 
and sin 2 26> E776 = 3.0 x 1CT 3 . The notation is the same as that in Fig. 1. 
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Figure 4: 

Predicted upper bound of CP violation in matter versus s 12 with = 163°. The 
notation is the same as that in Fig. 2. 
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